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' Abstract 

Q-f We have developed a framework for the Monte-Carlo simulation of the X-Ray Telescopes (XRT) and 

the X-ray Imaging Spectrometers (XIS) onboard Suzaku, mainly for the scientific analysis of spatially and 
spectroscopically complex celestial sources. A photon-by-photon instrumental simulator is built on the 
ANL platform, which has been successfully used in ASCA data analysis. The simulator has a modular 
structure, in which the XRT simulation is based on a ray-tracing library, while the XIS simulation uti- 

J> \ lizes a spectral "Redistribution Matrix File" (RMF), generated separately by other tools. Instrumental 

characteristics and calibration results, e.g., XRT geometry, reflectivity, mutual alignments, thermal shield 
transmission, build-up of the contamination on the XIS optical blocking filters (OBF), are incorporated 
as completely as possible. Most of this information is available in the form of the FITS (Flexible Image 
Transport System) files in the standard calibration database (CALDB). This simulator can also be utilized 
to generate an "Ancillary Response File" (ARF), which describes the XRT response and the amount of 
OBF contamination. The ARF is dependent on the spatial distribution of the celestial target and the 
photon accumulation region on the detector, as well as observing conditions such as the observation date 
and satellite attitude. We describe principles of the simulator and the ARF generator, and demonstrate 
their performance in comparison with in-flight data. 

Key words: Instrumentation: detectors — Telescopes — X-rays: general — Methods: data analysis 



1. Introduction 

A Monte-Carlo simulator is useful in characterizing a 
detector, and can relatively easily take into account many 
of the parameters which affect observations. Since the ul- 
timate goal of X-ray data analysis is to estimate the true 
time, energy and position of the incoming X-ray photons, 
it is quite important to predict precisely how the photons 
interact with the telescope and detector. A good simula- 
tor is therefore strongly required not only for instrumental 



calibration and proposal planning, but also for scientific 
analysis. Chandra and XMM-Newton also have good sim- 
ulators, named MARX 1 and SciSim 2 , respectively. 

The X-ray observatory Suzaku (formerly known as 
Astro-E2) is the fifth Japanese X-ray astronomy satel- 
lite (Mitsuda et al. 2006). It has been developed under a 
Japan-US international collaboration, and was launched 

1 http://space.mit.odu/ASC/MARX/ 

2 http:/ /xmm. vilspa.esa.es/scisim/ 
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on 2005 July 10. Five X-Ray telescopes are present, sen- 
sitive to soft X-rays below ~ 10 kcV (XRTs; Scrlcmitsos 
ct al. 2006). At the foci of four of the XRTs (XRT-I) 
are charge-coupled devices (CCD), known as the X-ray 
Imaging Spectrometers (XIS; Koyama et al. 2006); one 
(XRT-S) is combined with an X-ray calorimater known as 
the X-Ray Spectrometer (XRS; Kelley et al. 2006; XRS 
quit operation ~ 1 month after the launch). 

The Suzaku XRT is characterized by large collective 
area and relatively short focal lengths, compared with 
those of Chandra and XMM-Newton. In combination with 
these features, the low-earth orbit of Suzaku, where the 
particle background is low and stable, makes the non-X- 
ray background (NXB) of Suzaku much lower compared 
with Chandra and XMM-Newton. In addition, the XIS 
achieves good spectral resolution, especially at the low en- 
ergy range below ~ 1 keV with the backside-illuminated 
(BI) CCD for XIS1. The front-illuminated (FI) CCDs 
for XIS0, XIS2, and XIS3 exhibit about half of the NXB 
rate than XIS1 (and less at energies > 8 keV), so they are 
complementary. Therefore, Suzaku has a unique advan- 
tage for spectroscopic observations of spatially extended 
sources (Mitsuda et al. 2006). 

To achieve large collective area within the tight weight 
budget (1706 kg), the Suzaku XRT adopts the conical 
approximation of Wolter type I optics with 175 layers of 
the thin-foil-nestcd reflectors per quadrant (Serlemitsos et 
al. 2006). In return for the high throughput, it provides 
a moderate imaging capability of 2' half power diameter 
with a complex point spread function (PSF), as well as 
the energy-dependent vignetting effects common to X-ray 
telescopes. In addition, there exists spatially-dependent 
contamination on the optical blocking filters (OBF) of the 
XIS (Koyama et al. 2006). These XRT and XIS character- 
istics often make extended source analysis complicated, so 
it is crucial to prepare a tool in order to precisely evaluate 
the effect of complex telescope and detector responses. 

We developed a Monte-Carlo simulator of the Suzaku 
XRT/XIS system, which is incorporated into two practi- 
cal tools, the XIS simulator "xissim" and the "Ancillary 
Response File" (ARF) generator "xissimarfgen" . The sim- 
ulator is constructed on the "ANL" platform (§ 2.1), which 
is used for almost all of the processing and analysis soft- 
ware of Suzaku. 

While these tasks provide vast flexibility to the Suzaku 
XIS users, it is rather difficult to utilize them efficiently 
and appropriately. For example, there are more than 90 
parameters for both xissim and xissimarfgen. There are 
several issues and limitations that one should be aware 
of in running these tasks. This paper is aimed to clarify 
these things by explaining principles of the software and 
by demonstrating performance with practical examples. 
We have also tried to separate the 'calibration issues', 
which can be changed (usually improved) by calibration 
updates, from those originated in the design of the soft- 
ware itself. The quality of the calibration is out of scope 
for this paper, although some aspects are discussed briefly 
in §6. 

This paper is organized as follows. In § 2, we briefly 



show the strategy of the Suzaku software development, fo- 
cusing on the ANL platform and simulators. In § 3 and § 4, 
we describe principles of xissim and xissimarfgen, respec- 
tively. In §5, several notes on these tasks are described. 
In § 6, we demonstrate these tasks with three distinct ex- 
amples, the Crab Nebula, the North Ecliptic Pole (NEP) 
field, and Abell 1060. Finally, a summary is given in § 7. 
We also added three appendices which describe the coor- 
dinates definition, structures, parameters, and the output 
file formats, in detail. 

2. Software Development for Suzaku 

2.1. The ANL Platform 

When ASTRO-E software development started in 1995, 
the goal was a common software framework/platform 
which is used by realtime quick-look, data processing, and 
scientific analysis, both during pre-launch phase and af- 
ter the launch. To that end, it was necessary to provide 
a common programming environment where instrument 
team members can easily develop, maintain and updates 
softwares that they need. This framework/platform also 
must allow end-users to share these softwares. The frame- 
work must be easy to learn for instrument team members, 
who, spending most of the time in calibrating the instru- 
ments, do not necessarily have extensive programming ex- 
perience. Also, from the end users' point of view, it is de- 
sirable that those software tools developed based on this 
framework are maximally flexible and have an Ftools- 
like simple interface which is familiar to most X-ray as- 
tronomers. 

A software platform called "ASCA_ANL" , which had 
been developed for the ASCA satellite (Tanaka et 
al. 1994), fulfills these requirements. The ASCA.ANL 
platform mandates modular design of the analysis soft- 
ware to be built upon it, and makes the software prod- 
ucts easily configurable and reusable in components, so 
that software developers and end-users can share the same 
components for different purposes. This feature not only 
reduces code duplication, but also helps to quickly mature 
and refine the software. 

Indeed, ASCA_ANL fostered many practical tools in- 
cluding the instrument simulator SimASCA and the 
response generator SimARF. The advantages of the 
ASCA_ANL platform are demonstrated by original sci- 
entific research which would have been difficult without 
SimASCA and SimARF; e.g., spatial-spectral analysis of 
clusters of galaxies (Ikebe et al. 1996; Honda et al. 1996), 
systematic analysis of large volumes of X-ray surveys and 
the cosmic X-ray background (Ueda et al. 1998; Ueda 
et al. 1999; Kushino et al. 2002). The SimASCA and 
SimARF were very helpful to realize a specialized anal- 
ysis method in the analysis of spatially extended sources 
(which is not supported by standard analysis software), 
and to accurately compute complicated instrument re- 
sponses. 

On the other hand, however, ASCA_ANL and other rel- 
evant software were based on the functions and libraries 
used for realtime quick-look software which had been dc- 
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veloped by the instrument teams, independently from the 
official ASCA analysis software (Ftools). This resulted 
in two independent streams to calculate basic physical val- 
ues from the raw data, such as the pulse height corrected 
for the detector gain changes (known as "pulse invariant" , 
or PI, corresponding to the detected photon energy), and 
the sky and detector coordinates of events, which caused 
confusion in the scientific analysis of the ASCA data. 

Based on the ASCA experience, we adopted the "ANL" , 
i.e. a generalized version of the ASCA_ANL, as the soft- 
ware development platform for ASTRO-E and Suzaku. 
A brief history and concept of the ANL are described 
in Ozaki et al. (2006). At the same time, we devel- 
oped a mechanism to convert ANL software directly to 
Ftools, to ensure that the ANL tools used for calibration 
by instrument teams are equivalent to Ftools used for 
pipeline processing and scientific data analysis. Common 
FITS-read and -write ANL modules and functions were 
also developed to handle photon event files and the cal- 
ibration files in FITS format. 3 Now, almost all of the 
Ftools for Suzaku including xissim and xissimarfgen, re- 
leased from the Guest Observer Facility at NASA/GSFC, 
are developed in the ANL framework. 

2.2. History of Suzaku Simulators 

The development of the Suzaku simulator had started 
before the failure launch of the ASTRO-E on 10 Feb 
2000, especially for the bilinear 16 x 2 pixel XRS detec- 
tor (Kelley et al. 1999). The detector size was compara- 
ble with the angular resolution of the XRT (Kunicda et 
al. 2001; Shibata et al. 2001), and so the XRS PSF was 
undcrsampled. The energy resolution of the XRS was also 
very dependent on the count rate of each calorimeter pixel. 
Therefore, the XRT/XRS system simulator, xrssim, was 
required to estimate the flux coming to each XRS pixel, 
which was critical for proposal planning. This was also the 
case for the Suzaku XRS. The xissim task subsequently 
developed by replacing the XRS component of the simu- 
lator with the XIS component. 

The XRT ray-tracing part of the simulator has been sig- 
nificantly updated from the ASCA era. The code had been 
rewritten, by R. L. Fink (NASA/GSFC), from Fortran 
into C++, and the structure had been re-designed to 
utilize the mirror geometry and reflectivity files as sep- 
arate calibration FITS files. This ray-tracing code is now 
supplied as the "xrrt" library by the XRT team. It has 
been utilized for the performance improvement of the XRT 
(Misaki et al. 2005) and the design of the pre-collimator 
to suppress the stray-light (Mori et al. 2005). 

At present, xissim and xissimarfgen are publicly re- 
leased in the Suzaku Ftools and all the source code 
is available including the ANL itself and the xrrt li- 
brary. The latest version of the xissim package is 2006- 
08-26, which will be included in the next official release of 
the Suzaku Ftools for version 2.0 processing of Suzaku 
archival data scheduled in late 2006. All the calibration 

3 The calibration FITS files are released to the public from the 
NASA/GSFC guest observer facility, as part of the official cali- 
bration database (CALDB; George et al. 1991). 



information currently available is taken into account via 
the CALDB calibration database. The mkphlist, xissim, 
xissimarfgen, and xiscontamicalc tasks described in this 
paper are based on this version of the xissim package. 
The latest information on the xissim package is available 
at htt p : / / www-x. phys . met ro-u . ac . j p/ ~ ishisaki / xissim / . 

We also note that there is another Monte Carlo simu- 
lator for Suzaku, based on the Geant4 toolkit (Gcant4 
Collaboration et al. 2003) with ANL++ (Ozaki et 
al. 2006). This can simulate interactions of cosmic rays 
(both X/7-rays and particles) with the satellite materials, 
such as satellite structures, shielding around detectors, 
and the detectors themselves. The main purpose of this 
simulator is to study response of the Hard X-ray Detector 
(HXD; Takahashi et al. 2006; Kokubun et al. 2006), and 
the NXB models for both HXD and XIS. See Terada et 
al. (2005) and Ozaki et al. (2006) for details. 

3. Simulator: xissim 

The xissim task simulates the interaction of the incident 
X-ray photons with the XRT/XIS system, using the XRT 
ray-tracing library and a spectral "Redistribution Matrix 
File" (RMF; see §4.1) for the XIS, and generates a sim- 
ulated event file. The format of the generated event file 
is a stripped-down version of that created by the pipeline 
processing of a real observation, so that users can analyze 
the simulated data in the same manner as the real data. 
To perform the simulation, users need to take three steps. 
First, the spatial distribution and the energy spectrum 
of the celestial source to be simulated must be specified. 
Second, a list of incident photons from the source needs to 
be prepared as FITS file(s). An auxiliary tool, mkphlist, 
may be used for this purpose. Third, the photon FITS 
file is passed to xissim, which then performs a photon-by- 
photon simulation, and creates a file of events detected 
by the XIS. In the following subsections, we describe how 
xissim performs the simulation. 

3.1. Photon Generation 

An auxiliary task mkphlist generates a list of faked 
photons of an X-ray source from the model spectral 
energy and spatial distribution, photon flux (in pho- 
tons cm~ 2 s _1 ) in an arbitrary energy band, and the ge- 
ometrical area of the XRT (in cm 2 ), provided by users. 
A model spectral distribution file (which is specified by 
the qdp_spec_file parameter) must be in units of photon 
flux (photons cm~ 2 s _1 kcV -1 ) that can be easily pro- 
duced with standard software packages such as Xspec 
(Arnaud 1996). mkphlist requires celestial coordinates of 
the point source or a surface brightness map (FITS image) 
on the sky for the spatial distribution of the source. Either 
the number of photons or exposure time is needed to de- 
termine how many photons are to be generated. Users can 
also specify equal or random interval steps for the photon 
arrival time. The structure of mkphlist is explained in 
Appendix 2.1, and a list of major parameters and the for- 
mat of the photon file are summarized in table 5 and 8, 
respectively. Note that, by preparing an appropriate pho- 
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ton file, users can in principle simulate any source with 
any energy spectrum and/or any spatial distribution. 

3.2. Photon-by- Photon Simulation 

By taking into account the XRT and XIS response, xis- 
sim performs photon-by-photon simulation for given input 
photon file(s). It has the capability to read up to eight 
photon files simultaneously. Figure 1 shows the schematic 
structure of the simulation implemented in xissim. Since 
understanding the coordinate systems is essential, we in- 
clude the definitions in Appendix 1. 

First of all, the RA and DEC values in the photon file 
need to be converted to (9, </>), i.e. offset angle (') and 
azimuth angle (°), with respect to the XRT optical axis. 
This requires the satellite Euler angles {ea%, ea,2, ea^) (°), 
the observation date for an aberration correction (or par- 
allax correction, see Appendix 1), and the alignment pa- 
rameters in the telescope definition (teldef) file. Users 
can supply an attitude file (set of Euler angles as a func- 
tion of time) and a good time interval (GTI) file to take 
into account the wobbling of the spacecraft (See also § 5.2 
for the attitude wobbling). The photon_time column 
in the photon file usually starts from 0.0 s unless other- 
wise specified, and it is treated as the time offset relative 
to the GTI. Alternatively one may specify a fixed set 
of Euler angles and/or a fixed date. The aberration cor- 
rection can be disabled by setting the hidden parameter 
aberration=no (hidden parameters are not required when 
invoking an Ftools task). 

In the second stage, the geometrical area for a given 
photon is reduced by a factor of cos 8 due to the slanted 
incidence to the XRT. This factor is usually very close 
to unity, and had been neglected in the older version of 
xissim. This behavior can be controlled with the param- 
eter aperture_cosine, and is set to 'yes' by default in the 
present version. The photon flux is further reduced due 
to transmission through the thermal shield on the top of 
the XRT. Xissim then assigns a random location for each 
photon at the top surface of the XRT, where the pre- 
collimator is placed. The task traces the path of each 
photon inside the XRT (pre-collimator, primary and sec- 
ondary mirrors), using the XRT ray-tracing library, xrrt 
(Misaki et al. 2005; Mori et al. 2005), using the XRT geom- 
etry and reflectivity as described in the ray-tracing code 
and the calibration files. After the ray-tracing, some pho- 
tons may be absorbed and disappear, while others reach 
the focal plane. 

A fraction of the photons that have reached the focal 
plane are absorbed by the contamination on the OBF. 
The thickness of the contamination is time- and detector- 
position-dependent (Koyama et al. 2006), and their de- 
pendence is given by a calibration file supplied by the 
XIS team. Xissim computes the transmissivity at a given 
time and position using this calibration file. The position 
of the photon on the detector is again calculated by the 
alignment parameters in the teldef file. 

Finally, the simulated photons reach the detector (in- 
cluding both OBF and CCD), where the detection prob- 
ability is determined using the RMF of the XIS. The 
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Fig. 1. Schematic structure of the simulation. 

XIS RMF contains the transmission of the OBF and the 
quantum efficiency of the CCD, as well as the spectral re- 
distribution matrix from energy to PI. The line response 
function of the XIS CCDs is primarily a Gaussian distri- 
bution but it also includes other features such as escape 
ratios and tails that deviate from a Gaussian. Photons 
that have passed the test for detection are recorded as 
X-ray events, 4 and their PI values are determined from 
the incident photon energy by random choices accord- 
ing to the energy redistribution probability in the RMF. 
The Suzaku XIS detectors do not exhibit significant posi- 
tional dependence in the energy resolution after the CTI 
(charge transfer inefficiency) correction while the energy 
resolution is known to degrade with time. Users should 
supply an appropriate RMF corresponding to the obser- 
vation date, which can be generated by a separate task, 
xisrmfgen. 

Note that the current version of xissim does not consider 
the NXB, bad CCD columns, event pile-up, event grade, 
nor CCD exposure frames. Although the output event 

4 We shall call 'photon' during the simulation, which becomes 
'event' after the detection. 
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files contains the same major columns as the event files of 
the real data, the STATUS and GRADE columns are filled 
with 0, while the pha column has the same value as the 
pi column. 

The ANL module structure and parameters of xissim 
is explained in Appendix 2.2 in detail. If one has the 
ANL programing environment available, he/she may add 
his/her own modules to the simulator. It is also easy to 
replace a module, e.g. if a module is available that more 
precisely simulates the CCD detection process, then this 
module can be substituted for the SimASTEJCISRMFsim 
module which utilizes a ready-made RMF. This is one of 
the great benefits of the ANL. 

3.3. Calibration Files 

Table 1 summarizes the list of calibration files used by 
xissim. The file specified by the leapfile parameter is the 
leap second file, and is required to compute the mission 
time (or Suzaku time), defined as accumulative seconds 
since 2000 January 1, 00:00:00 (UTC). In fact, the de- 
fault value of the leapfile parameter is set to a special 
keyword of "caldb" and is a hidden parameter. By in- 
stalling CALDB and properly setting the environmental 
variables, the xissim task automatically searches the most 
recent calibration file for this category, i.e. 'Content Name' 
= leapsecs. This is also applicable to other parameters 
in table 1. 

The shieldfile, mirrorfile, reflectfile, and backproffile pa- 
rameters are used for the XRT simulation. There are four 
FITS extensions in the mirrorfile to describe the geome- 
try of each XRT, and three extensions are present in the 
reflectfile corresponding to materials of the reflection sur- 
face. As described in Appendix 1, teldef is used to describe 
the mutual alignments between XRT and XIS, as well as 
among the XIS sensors and the spacecraft. The contam- 
ifile describes the energy, time, and position dependence 
of the contamination on the XIS OBF. 

4. Ancillary Response Generator: xissimarfgen 

The xissimarfgen task generates a Suzaku XIS ARF 
based on user-defined conditions, such as an arbitrary 
shape of the X-ray emitting region and event extraction 
regions. Xissimarfgen does so by simulating photon de- 
tections at each energy. It then calculates the detection 
efficiency in a user-defined event accumulation region. 
Since it utilizes a Monte-Carlo simulation, users need to 
simulate enough photons to avoid counting statistics er- 
rors. It can refer to the attitude file to reflect the change 
of effective area due to the attitude wobbling. The final 
ARF is in the standard FITS format, so that users can 
use Xspec or other standard fitting packages for spectral 
analysis. 

4-1. Principle of ARF Calculation and Limitations 

The ARF is utilized for spectral fitting combined with 
an RMF. See George et al. (1992) for detailed format 

5 We shall use 'accumulate' for the simulated events, and 'extract' 
for the real observed events. 



of these files. The RMF is represented by an (m x n) 
matrix R(Ei,PIj), where E (keV) denotes the energy 
and PI (channel; hereafter chan) denotes the pulse in- 
variant, with l<i<m and 1 < j < n. Regarding the XIS, 
m = 7900, n = 4096, E l = 0.201 keV, E m = 15.999 kcV, 
PI i = chan, and PI n = 4095 chan for the nominal RMF. 
The ARF is represented by an m-dimensional vector which 
we denote as S A(Ei) (cm 2 ), where S = 1152.41 cm 2 rep- 
resents the geometrical area of the XRT. The goal of the 
spectral fitting is to find a model spectrum, M(Ei) (pho- 
tons cm~ 2 s _1 kcV -1 ), which fits the observed spectrum, 
V(PIj) (count cm~ 2 s _1 ). The response and model spec- 
trum are convolved, i.e., 

m 

M(PIj) =S^ AE i A (Ei) R{Ei,PIj) M(Ei), (1) 

where AE (keV) is the energy bin width, and A4(PIj) 
and T>(PI j) are compared. As one can see easily from 
this formula, [ A(Ei) R(Ei, PI j) ] represents an expected 
spectrum for the monochromatic X-ray of E = Ei keV, and 

n 

A(Ei) ~y^R(Ei,PIj) represents the detection efficiency at 
E = E, keV. 

Thus, calculating the ARF is reduced to the computa- 
tion of the detection efficiency at each energy step, £j, 
of the RMF, a job well-suited for a Monte-Carlo simu- 
lation. For a given input N[ n counts of monochromatic 
X-ray photons at E — Ei keV, the simulator predicts A^et 
detected events and then the detection efficiency is simply 
A(Ei) = Adet/Am- However, one must be very careful be- 
cause the detection efficiency, namely N^ct, is influenced 
by many factors: first of all, the accumulation region of 
the event on the detector, and the spatial distribution of 
the celestial sources assumed on the sky. It is also affected 
by the observational conditions, such as the satellite Euler 
angles, the date of the observation due to the thickness of 
the XIS contamination and the parallax correction, etc. 
The quality of the calibration and/or the Poisson statistics 
can also impact Adct- It is therefore important that one 
must reproduce the user-selection and the observational 
conditions of the real data as much as possible in the sim- 
ulation. One must also take care to perform a simulation 
such that the photon statistics are sufficiently better than 
the statistics of the real observation. 

In fact, the spatial distribution on the sky is sometimes 
complex and/or extended on a scale larger than the tele- 
scope FOV. Thus the accuracy of the spatial model can 
become a major cause of systematic error in the estima- 
tion of the detection efficiency, which leads to uncertainty 
in the source flux. For example, if one assumes a more 
core-concentrated image than in reality, more photons will 
be simulated to arrive at the detector, which will over- 
estimate the detection efficiency. One can test the as- 
sumed spatial distribution on the sky by comparing the 
real observation image and the simulated one. 

There is also another limitation due to the spectral 
fitting procedure itself. In the conventional spectral fit- 
ting package (e.g., Xspec vll or before), one can choose 
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Table 1. List of calibration files used by xissim 



Parameter File Name * 



Content Name T Description 



leapfile leapsec_010905 . f its 

shieldfile ae_xrta_shield_20061129 .fits 
mirrorfile ae_xrtiV_mirror_20060710 .fits 



reflectfile ae_xrta_reflect_20060710.f its 



backproffile ae_xrta_backprof .20060719 .fits 
teldef ae_xiiV_teldef .20060125 .fits 

contamifile ae_xiiV_contanii_20060525 .fits 



LEAPSECS 

FTRANS 

GEOMETRY 

GEOMETRY 

GEOMETRY 

GEOMETRY 

REFLECTIVITY 

REFLECTIVITY 

REFLECTIVITY 

BACKPROF 

TELDEF 

CONTAMI_GROWTH 
CONTAMI_TRANS 



Table of times at which leap seconds occurred 

XRT thermal shield transmission 

XRT mirror geometry 

XRT obstruction geometry 

XRT quadrant geometry 

XRT pre-collimator geometry 

XRT mirror foil front surface reflectivity 

XRT mirror foil back surface reflectivity 

XRT pre-collimator surface reflectivity 

XRT foil backside scattering profile 

Telescope definition file 

XIS OBF contamination growth curve 

Template transmission vs energy for the contaminant 



* N represents 0, 1, 2, or 3 respective to the XIS sensor. 

t The CCNMmmn keyword in the FITS header, and the CAL.CNAM column in the CALDB index file. 



only a single response matrix (ARF + RMF) for an ob- 
served spectrum in the spectral fitting. 6 For example, an 
observed spectrum may contain thermal emission which 
obeys an oval surface brightness profile, as well as the 
cosmic X-ray background (CXB) spectrum of a T ~ f .4 
power-law which extends nearly uniformly on the sky. The 
ARF response for the oval surface brightness is different 
from that for the uniform-sky emission, hence one can- 
not fit the observed spectrum with the thermal model + 
power-law model in a usual way. Strictly speaking, the 
energy spectrum should be the same at every point in the 
assumed spatial distribution on the sky in order to con- 
duct spectral fitting with a single ARF + RMF response. 

4-2. Implementation of ARF Calculation 

As described in §3.2 and shown in figure 1, the XIS 
RMF takes care of the OBF transmission and the quan- 
tum efficiency of the CCD, hence the XIS ARF should 
consider other factors for the detection efficiency, namely, 
the thermal shield transmission, XRT effective area, trans- 
mission of the OBF contaminant, etc. Detailed explana- 
tion of structure, parameters, and the output ARF format 
are given in Appendix 2 and 3 

It reads a number of parameters which specify the 
simulation conditions (table 7), and (1) determines en- 
ergy steps to calculate detection efficiency; (2) generates 
monochromatic photons (or quasi-monochromatic within 
the narrow energy range) until the user-specified condition 
on the photon statistics is fulfilled at each energy step; (3) 
conducts the ray-tracing simulation for each photon; (4) 
counts up the number of detected events at each energy; 
(5) records the detection efficiency at each RMF energy 
bin to the output ARF(s) by interpolating the simulation 
result; (6) continues to the next energy step and loops to 
step (2). 

Note that the energy step determined in step (1) is usu- 
ally not same as the RMF energy bin, because the com- 

6 This restriction no longer holds in the latest release of Xspec 
vl2, which allows different model components to have their own 
response. 



putation time would be very long to conduct photon-by- 
photon simulations in standard XIS RMF 2 eV steps up 
to 16 keV. Interpolation is therefore required in step (5). 
In addition, the XRT effective area usually changes only 
gradually with energy except for several characteristic en- 
ergies such as the Au-M, Au-L, and Al-K edges, 7 so that 
we may often choose sparse energy steps. This feature can 
save the computation time effectively. 

In the calculation of the detection efficiency, the 
three major factors of (i) transmission of the XRT 
thermal shield, (ii) effective area (cm 2 ) of the XRT, 
and (iii) transmission of the XIS OBF contaminant 
are treated separately. They are also written in sepa- 
rate columns in the resultant ARF as SHIELD_TRANSMIS, 
XRT_effarea, and CONTAMI_transmis (table 10). The 
resultant detection efficiency times the geometrical area, 
S A(Ei) (cm 2 ), is written in the SPECRESP column, i.e., 

SPECRESP = SHIELD_TRANSMIS X XRT_EFFAREA X CON- 

tami_transmis. Note that (i) and (iii) arc supplied in 
the calibration files (table 1) in fine energy steps of ~ eV, 
whereas (ii) is usually calculated in more sparse energy 
step. By separating these factors, one can obtain a good 
quality ARF even in a sparse energy step for the simula- 
tions, and moreover, one may remove, scale, or multiply 
the CONTAMI.TRANSMIS factor afterwards. The xiscon- 
tamicalc task is provided to do this kind of the ARF ma- 
nipulation. 

Note that the thickness of the OBF contaminant is po- 
sitionally dependent. It is therefore required to know 
the spatial distribution of photons 8 falling on the OBF 
at each RMF energy bin in order to evaluate the CON- 
TAMI.TRANSMIS factor. This energy dependence of the 
photon distribution is also determined by interpolation, 
which incurs additional calculation time when the simula- 

7 The front surface of the XRT reflector is coated with gold and 
its substrate is made of aluminum. The pre-collimator is made 
of aluminum, too, so that the Al-K edge appears in the large- 
offset-angle response of the XRT. 

8 This distribution is approximated in xissimarfgen by a DET 
coordinate image binned prior to applying absorption due the 
XIS contamination. 
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tion energy step is much wider than the RMF energy bin. 
It is not easy to estimate the true photon distribution from 
the real observation data, because the observed image is 
affected by the XRT vignetting and the OBF contami- 
nant, both of which are energy dependent. Vignetting is 
a more severe effect in the higher energy band, and the 
OBF contamination is severe in the lower energy band. 
In addition one must subtract background to utilize the 
observed image. The combined energy and spatial depen- 
dence of the XIS contamination is considered in the ARF 
generator rather than the RMF generator, for this reason. 

At each simulation energy, in fact, the A(E) value is 
calculated using the weighted sum of events, iV w , instead 
of A dot , as, 

N in (E) 

A(E) = N W (E) I Ni n (E) = ME)/N in (E), (2) 

k=l 

in which Wk{E) denotes the WEIGHT value (see 
Appendix 2.2) of each simulated photon at the energy 
of E keV. As mentioned above, the resultant A{Ei) val- 
ues at the RMF energy bin are calculated by interpola- 
tion, complicated somewhat by when the transmission 
of the OBF contaminant is considered. Here, we de- 
fine I = iNDEXi (see Appendix 3.3 for index*). There 
are Ni n (El) photons with weight without contamination 
represented by Wk(E[) and energy a little below Ei, and 
A; n ) photons with w/-'(E' l+1 ), energy a little above 
Ei, i.e., E[<Ei< E' l+l . The transmission of the OBF con- 
taminant is calculated for each of the simulated photons, 
as t k {Ei , PHOTON_TiMEfe , DETXj, , DETYfe ) . Note that the 
energy of each simulated photon, E[ = PHOTON_ENERGYfc, 
has been replaced by the energy of the RMF bin, Ei. 
Thus xissimarfgen computes the final detection efficiency, 
A(Ei), with contamination, as 



A(Ei) = Si 



Tk w k 



E 



(3) 



by an interpolation. The definitions of Si and T* are given 
in eqs. (A2) and (A3). 

It also calculates the relative error of A(E) at each sim- 
ulation energy, and the interpolated values are stored in 
the relerr column of the output ARF (table 10). This 
column is useful to judge the photon statistics is sufficient 
for the ARF calculation. The relative error is calculated 
as, 



R.ELER.R, — 



A in - A dct 
Ai n A d ct 



ct 



1 



(4) 



if A in A dct (A in - A dot ) 7^ 0, otherwise relerr = 1.0. The 
derivation of this formula is a little tricky, because we 
know the detected count A de t and the undetected count 
Ai„ — A dot in the simulation, and both are considered 
to follow the Poisson statistics. Since A(E) is expressed 
as A(E) = A dot /A in = 1 - (A in - A dot )/A in , the error of 
A(E) can be evaluated in two ways, 8A\ = VAdct/Am 
or SA2 = VAi n — Ad c t/Ai n . We therefore defines the 

-(6A 2 )-yA = 



V(N ia - A dct )/A in /A dct - eq. (4). 
5. Notes 

In this section, we describe several notes on xissim and 
xissimarfgen. § 5.1 applies to both tasks, and others apply 
mainly to xissimarfgen. 

5.1. Notes on Random Numbers 

The quality of the random number generator to be used 
can affect the quality of the Monte-Carlo simulation re- 
sults. A good random number generator should include 
a very long cycle, fast computation, and wide significant 
bits, xissim and xissimarfgen use an internal random num- 
ber generator in the astetool library, 9 utilized by all mod- 
ules. This generates double precision floating point values 
in the range of < r < 1.0 based on the Tausworthe method 
(Tausworthe 1965). The generated random number has 62 
significant bits (~ 4.6 x 10 18 ) and its cycle is estimated to 
be about 2 250 ~ 10 75 . These parameters are significantly 
wider and longer than the usual random number function, 
int rand(void), implemented in the standard C library. 

Its code is machine independent, and it reproduces ex- 
actly the same series of random numbers as long as the 
rand_seed and rancLskip parameters are the same. It is rec- 
ommended to set a prime number (except 2) to rancLseed 
for good randomization. The default value of rand_seed for 
the simulation tasks is 7. They record the number of ran- 
dom numbers generated in the simulation to the output 
event file, as the RANDNGEN keyword in the FITS header. 
One may re-continue the simulation with the same series 
of random numbers by setting the rand_skip parameter to 
its value. However this code is not multi-thread compli- 
ant, which may need to be upgraded in the future for 
faster (i.e., distributed) simulations. 

5.2. Notes on Accumulation Region 

There is a difference between specifying the accumula- 
tion region in SKY coordinates versus DET coordinates. 
This may be ignored only when the attitude wobbling and 
the parallax correction are negligible. The accumulation 
region is fixed on the CCD when it is specified in DET co- 
ordinates. On the other hand, it moves around the CCD 
when specified in SKY coordinates, according to the atti- 
tude wobbling. In both cases, the celestial target moves 
around the CCD and is affected by the vignetting effect of 
the XRT, also due to the attitude wobbling. Xissimarfgen 
can treat both situations correctly, as far as the supplied 
attitude file is reliable, so that one should select the re- 
gion_mode parameter to match the extraction method of 
the real observation spectrum. 

It is known that there is an unexpected attitude wob- 
bling of ~ 0.5' due to thermal distortion of supporting 
structure (Serlemitsos et al. 2006), however this effect is 
not included in the present attitude file. This situation 
will be improved in near future by a dedicated Ftool, 
the aeattcor task. Until then, it is recommended to avoid 



relative error as 5 A/ A = 1/^(SA 1 ) 



See also Appendix 1 for the astetool library. 
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using too small of an accumulation radius (r < 3'). One 
may check this effect by changing the accumulation radius 
and test whether the fit results are affected significantly. 
Alternatively, one may track the position of the PSF core 
on the CCD for bright point-like source targets. 

It is also notable that the background files for the XIS 
currently released, which are a collection of events when 
the XRT was pointed to the night (non-sunlit) Earth, do 
not support event extraction in SKY coordinates. This 
situation will be improved in the future. One may extract 
the background from the outer ring of the target, how- 
ever, this region also contains the outskirt of the PSF of 
the main target, CXB, and the instrumental background, 
which have a small dependence on the detector position. 
The former two effects can be evaluated by xissimarfgen, 
and the latter can be tested with the released background 
file. 

5.3. Notes on Flux Normalization 

Because the detection efficiency defined in cq. (2) is 
considered for all the input photons coming from every- 
where in the supplied source image, the normalization of 
the flux in spectral fitting gives the value integrated over 
the whole region of the source image. Therefore, if one 
generates a uniform-sky ARF with source.rmax = 20' to 
fit the CXB spectrum, the fit gives the flux from the 
7r • source_rmax 2 = 1257 arcmin 2 sky area, then the user 
needs to divide the flux by this area to convert it to a 
surface brightness. 

Other cases can similarly be complex, e.g., an analysis 
of a cluster of galaxies. Extracting spectra from annular 
rings centered on the cluster core is frequently performed 
in the cluster analysis. Here, we assume that the clus- 
ter emission spectrum is identical everywhere on the sky, 
and only the normalization of the flux decreases as the 
distance from the cluster core increases. We also assume 
that the spatial distribution of the cluster on the sky can 
be perfectly predicted, which has been supplied to xissi- 
marfgen as the source image. Then the fit results for each 
ring should give the same flux, while the observed count 
per unit area decreases as the ring radius increases, since 
the flux for the whole cluster is calculated for each fit. 

If one gets different fluxes for each fit, then this is the 
1st order approximation of the correction factor to the as- 
sumed source image at each ring. It is often desired to 
derive the flux only coming from each ring. To help with 
this kind of task, there is a keyword, SOURCE_ratio_REG, 
written in the output ARF (table 11). This keyword holds 
the ratio of the source image inside the specified accumula- 
tion region for the ARF, which has been calculated during 
the simulation. By multiplying this factor by the obtained 
flux, the user can calculate the flux in that ring. 

5.4- Notes on Computation Time and Memory 

The code of xissimarfgen is designed to conduct the 
computation of an ARF as efficiently as possible in both 
time and memory, although it still requires a significant 
amount of both. The simulation code has been tuned 
for speed; it reads all the required information including 



Table 2. Examples of computation time. 



CPU Time * Ap Simulation Condition * 

(A) 179.5 s — Full, source_mode=SKYFlTS 

(B) 109.9 s -69.6 s -s- contamifile=NONE 

(C) 101.0 s -8.9 s -> aberration=no 

(D) 92.3 s -8.7 s -> sourcejnode- j2000 



* Total number of CPU-seconds that the process used directly in 
user mode, measured by the Linux time command, 
t Difference of time compared with the preceding line, 
t Parameters of num_photon=100000 and estepfile=SPARSE (55 en- 
ergy steps) are common to all. 

the attitude into memory before the simulation. Searches 
of tables such as reflectivity, transmission, spatial and 
spectral distributions are accelerated by adding an in- 
dex. Several functions cache previous values to skip re- 
dundant calculations especially when the photon energy 
and/or time is similar to the previous ones. In addition, 
the binary distribution of the xissim/ xissimarfgen pack- 
age is compiled with fast C compilers using the highest 
optimization option. 

The required memory is usually around 130 MB, hence 
recent machines can easily run xissimarfgen task in mem- 
ory. The actual calculation time is very dependent on the 
simulation energy step, the photon statistics, as well as the 
computer platforms. Table 2 shows examples of compu- 
tation time on the AMD Athlon™64 2.4 GHz CPU with 
a 64-bit Linux OS. The example (A) is the ARF shown 
in § 6.1, in which full observational features are taken into 
account, and a Chandra image (1800 x 1800 pixels) was 
supplied for the sourceJmage with source_mode=SKYFiTS. 
Parameters of num_photon=100000 and estepfile=SPARSE 
(55 energy steps) were chosen, so that 5.5 x 10 6 photons 
were simulated for the ARF calculation. 

The time difference between (A) and (B) indicates that 
significant fraction of time (~ 70 s) was consumed in the 
calculation of the XIS contamination. However, this time 
is only proportional to m x Adot and does not depend on 
the simulation energy step, hence it should be acceptable. 
The parallax (aberration) correction also needs the non- 
negligible cost of ~ 9 s, which is proportional to N[ n = 
num_photon. Similar time is needed for the randomization 
in the spatial distribution of the Crab nebula, as seen in 
(C) — > (D). The consumed time in the XRTsim module 
is also displayed by the ANL, and it was 75.9 s for (D). 
This indicates that the ray-tracing code can perform the 
simulation of a single X-ray photon in less than 15 /is on 
this machine. 

6. Demonstration 

In this section, we demonstrate how xissim and xissi- 
marfgen work using three distinct examples: the Crab 
nebula as a calibration source and a quasi-point-like 
source in §6.1, the NEP field as a "blank sky" in §6.2, 
and Abell 1060 as an spatially extended source in §6.3. 
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6.1. Crab Nebula 

First, we present the case of simulating the Crab neb- 
ula which is the main X-ray calibration source for effec- 
tive area calibration (Toor & Seward 1974; Seward 1992; 
Kirsch ct al. 2005). Small in angular scale, it has a 
complex spatial structure as seen in figure 2 (a) of the 
Chandra image (Weisskopf et al. 2000). With respect to 
the surface brightness map, we adopted this image, be- 
cause Chandra's X-ray telescope, HRMA, has much supe- 
rior angular resolution of ~ 0.5". We further compensated 
it manually for a point-like emission from the neutron star 
(K. Mori priv. comm.). 

Wc made a photon list by supplying the image to 
mkphlist, and ran xissim with it. The simulated Suzaku 
image of the Crab nebula is shown in figure 2 (c) . For com- 
parison, we also present the simulated image for a point- 
like source in figure 2 (b). The simulated Crab image ap- 
pears as a smoothed PSF with the extent of the complex 
surface brightness profile of the Crab nebula. Figure 2 (d) 
shows the real observation image taken with the Suzaku 
XIS0 detector. The global extent of the Crab image is con- 
sistent with that of the simulated image. The anisotropy 
in the azimuth direction in the real image is due mainly 
to the complex PSF shape of the actual XRT, which will 
be more accurately reproduced by future improvement of 
the mirror geometry file (mirrorfile). Once the calibration 
file is updated, xissim can reflect it automatically via the 
CALDB. A narrow groove crossing the central area from 
east to west is due to a bad CCD column. The out-of- 
time events, which broadly spread on both the east and 
west sides, are also seen along the direction of the signal 
transfer from imaging area to frame-store region. These 
features are not implemented in the current version of xis- 
sim. 

Using the Chandra image, we also generated an ARF 
for XIS0, and it is plotted in figures 3 (a) and (b). 
The SPECRESP (black) and XRT_effarea (green) columns 
are plotted in figure 3 (a), and the CONTAMI_transmis 
(black) and SHIELD_transmis (green) columns are plot- 
ted in figure 3 (b). The 90% confidence range of the 
SPECRESP is also drawn by cyan lines in figure 3 (a) . Full 
observational conditions, namely attitude, gtifile, aberra- 
tion, and contamifile, are considered in the ARF generation 
with num_photon=100000 and estepfile=SPARSE (55 en- 
ergy steps). The accumulation radius is 6 mm — 250 pixel 
~ 4.34' in the SKY coordinate. 

For comparison, we plot the nominal ARF in CALDB 
without contamination, ae_xi0_xisnom6_20060615 . arf , 
in red line. In fact, the nominal ARF was also generated 
by an older version of xissimarfgen, and the calibration 
files were not changed between the two versions. However, 
the nominal ARF is calculated with much denser energy 
step (2 eV steps below 4 kcV, and at most 10 eV steps 
above 4 keV, with 3450 energy steps), and 4 times higher 
photon statistics (num_photon=400000). Although slight 
jerks are seen in black and green lines in figure 3 (a), 
these two ARFs are quite consistent. Discrepancy in the 
lower energy range is due to the XIS contamination, which 



is plotted by a black line in figure 3 (b). Therefore, in 
the spectral fitting, the Crab nebula can be treated as a 
point-like source with Suzaku XIS, if the extraction radius 
is large enough (r ~ 6 mm) and the spacecraft attitude is 
stable. We also note that the nominal ARFs give flux con- 
sistent with that obtained by Toor & Seward (1974) within 
<~ 2% for all the XIS sensors (Serlemitsos et al. 2006). 

6.2. NEP Field 

The NEP field is an archetypal "blank field" , where no 
X-ray bright objects exist. In such a region, the X-ray 
background, including both the extra-galactic (Brandt & 
Hasingcr 2005) and the Galactic components (Snowden 
et al. 1995), is the dominant X-ray source. The X-ray 
background can be treated as almost uniform distribu- 
tion, hence we tested the uniform-sky ARF generated by 
xissimarfgen in this field. 

We created an ARF assuming a uniform distribution 
for the source from a circular region with a radius of 20' 
(see caption of figure 4 for details of the parameters) , and 
fitted the observed spectrum with it. In figure 4 (a), the 
effective area, SPECRESP, of the obtained ARF is displayed 
in comparison with that for a point source. One can see 
that the effective area is relatively smaller in the higher 
energy band (> 7 keV) , due to the vignetting effect of the 
XRT. After subtracting the NXB contribution estimated 
using the night Earth database (§5.2), the spectrum can 
be well fitted with a power-law model representing the 
CXB and one or two thin-thermal plasma models repre- 
senting local Galactic thermal components, as shown in 
figure 4 (b) . The photon index and the surface brightness 
of the power-law component are consistent with the pa- 
rameters reported so far (Gendreau et al. 1995; Kushino 
et al. 2002). See Fujimoto et al. (2006) for the details of 
the analysis. This result, demonstrates that xissimarfgen 
properly generates the ARF for the uniform-sky emission. 

6.3. Abell 1060 Cluster of Galaxies 

Finally, we present an example of the Abell 1060 cluster 
of galaxies observed with Suzaku. Scientific results will be 
published by K. Sato et al. in preparation. Abell 1060 is a 
circular and nearly isothermal (~ 3 keV) cluster of galaxies 
(Tamura et al. 2000; Furusho et al. 2001; Hayakawa et 
al. 2004; Hayakawa et al. 2006) and is suitable for testing 
the ARF for extended sources. There were two pointings 
performed with Suzaku at the central region and the ~ 
20' east offset region, as shown in figure 5 (a). These 
observations were conducted at the end of November 2005, 
when the XIS contamination was already significant and 
was starting to saturate. 

The observed spectrum is assumed to contain (a) thin 
thermal plasma emission from the intra cluster medium 
(ICM), (b) local Galactic emission, (c) CXB, and (d) 
NXB. We can estimate (d) using the night Earth 
database mentioned in §5.2, and can subtract it from the 
observed spectrum. As demonstrated in § 6.2, the spec- 
trum of (b) can be represented by the (apec + apec) model 
with I solar abundance, and that of (c) has a shape of 
absorbed power-law with r ~ 1.4. However, we cannot 
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Fig. 2. Observed and simulated images of the Crab nebula, (a) Chandra image (Weisskopf et al. 2000), corrected for pile-up. (b) 
Simulated Suzaku image for a point-like source, (c) Simulated Suzaku image using (a) as the spatial distribution, (d) Observed 
Suzaku XISO image smoothed with a Gaussian of <r = 4 pixel ~ 4". The image width of Chandra is 2.9', while those of Suzaku are 
11.6'. See text for details. 
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Fig. 3. Plots of ARF columns generated for the Crab nebula, with num_photon=100000 and estepfile=SPARSE (55 energy steps) 
using the Chandra image in figure 2 (a) as sourcejmage. The accumulation radius is 6 mm = 250 pixel ~ 4.34' in the SKY coordinate, 
(a) black: SPECRESP, green: XRT_EFFAREA, cyan: 90% confidence range of SPECRESP, namely, SPEORESP ± 1.65 X RESPERR, although 
it is almost hidden by the overlaid black line. Red line indicates the nominal ARF in CALDB without contamination, (b) green: 
SHIELD_TR ANSMIS , black: CONTAMI_TRANSMIS. 




Fig. 4. Plots of ARF columns generated for the NEP field, with num_photon=400000, estepfile=DENSE (2303 energy steps), 
source_mode=UNlFORM, and SOURCE_RMAX=20'. The accumulation region is all the XIS1 CCD including the calibration source, 
(a) black: SPECRESP, blue: SPECRESP / C0NTAMI_TRANSMIS, cyan: 90% confidence range of SPECRESP. Solid red line indicates the 
nominal ARF in CALDB without contamination, and dashed red line shows it multiplied by 0.1. (b) Black line in the upper panel 
represent the NEP field spectrum observed with Suzaku in the "stable" period (Fujimoto et al. 2006). The NXB is subtracted, and 
it is fitted by the [ apec (cyan) + apec (orange) + wabs X power-law (blue) ] model in XSPEC 11.3.2t indicated by green line. The 
estimated NXB spectrum is overlaid in red line. Fit residuals in units of a are shown in the lower panel. 
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fit the observed spectrum directly by a sum of (a) + (b) 
+ (c), because the spatial distribution of these three are 
different on the sky, as described in §4.1. 

We therefore adopted the following method for the spec- 
tral analysis of Abell 1060. We extracted several spectra 
from annular regions centered on the cluster core, and 
here we show two samples of the innermost region at the 
projected radius of 0-2' from the central observation, and 
the outermost region of 17-27' from the offset observation, 
as representatives. We generated two different ARFs for 
each spectrum, A v (Ei) and A B (Ei), which respectively 
assume the uniform-sky emission and the ICM surface 
brightness profile obeying an analytical model obtained 
with the XMM-Newton data. 

As described in §4.1, it is important that the assumed 
spatial distribution on the sky well agrees with the actual 
data in the calculation of the ARF response. We therefore 
compared the observed images with the simulated ones in 
figure 5 (b) and (c). The 1-4 keV energy range was chosen 
so that the distortion of the image due to the XIS contam- 
ination and the XRT vignetting was not severe. In this 
energy range, the Galactic component (b) is almost negli- 
gible, whereas the CXB and NXB components cannot be 
neglected especially in the offset observation. The NXB 
component (red line) is estimated from the night Earth 
database. A small ACTY (= DETY-1 for XISO, see fig- 
ure 8) dependence of the NXB intensity is seen, reflecting 
the dwell time at the frame-store region of the CCD. The 
CXB component (blue line) is estimated by the xissim 
simulation, assuming the uniform sky and the previous 
ASCA results of the CXB intensity (Kushino et al. 2002). 
The vignetting effect is seen in the CXB counts, hence 
the count rate slightly drops at the CCD rim. After the 
subtraction of the estimated CXB and NXB components, 
the observed distribution of the cluster (black crosses) is 
fairly well reproduced by the xissim simulation of the clus- 
ter emission (green line), although a small asymmetry is 
observed for the real cluster in the central observation. 

Figures 6 (a)-(d) show the latter kind of ARFs, A B (Ei), 
for both regions. Figures 6 (a) and (b) correspond to 
the extraction regions of 0-2' and 17-27', respectively, in 
the DET coordinate, and the calibration source area (top- 
left and bottom-right for XIS1) is also excluded in (b). 
Although the accumulation area is smaller for (a) than 
(b), the calculated effective area is much larger for (a) 
than (b) as seen in figure 6 (c), plotted in black and red 
lines, respectively, due to the assumed surface brightness 
profile. One can see the position dependence of the XIS 
contamination (thinner towards the CCD edge) is treated 
appropriately as seen in figure 6 (d). 

Denoting the spectra of (a), (b), (c), and (d) as 
M ICM (Ei), M OAL (Ei), M CXB (Ei), and M" XB (PIj), the ob- 
served spectrum can be expressed by a sum of, 

a b M icm + A v M G AL + A v M cxb + 

where the operator ® denotes the transformation defined 
by eq. (1). It is known that the CXB spectrum, M CXB , 
is fairly constant over the sky except for difference in 
the neutral hydrogen column density, Ah , for absorption, 



whereas the local Galactic emission, M GAL , may vary from 
field to field by more than an order of magnitude (Kushino 
et al. 2002). 

Considering this situation, we assumed a power-law 
spectrum for the CXB with the values by Kushino et 
al. (2002), r = 1.4 and S x = 5.97 x 10~ 8 erg cm" 2 s" 1 sr" 1 
(2-10 keV), 10 measured with the ASCA GIS (Ohashi et 
al. 1996; Makishima et al. 1996). The neutral hydro- 
gen column density was fixed to Ah = 4.9 x 10 20 cm~ 2 
(Dickey & Lockman 1990). We calculated the estimated 
contribution of the CXB, A v £g> M CXB , using the fake com- 
mand in Xspec. This contribution for each region is 
indicated by blue crosses in figures 7 (a) and (b) for 
XIS1. We subtracted the CXB contribution from the 
observed spectrum as well as the estimated NXB spec- 
trum. The XIS1 (red) and FI (XIS0+XIS2+XIS3; black) 
spectra in figures 7 (a) and (b) denote those after the 
CXB and NXB subtraction. We then fitted the XIS1 
and FI spectra simultaneously for the offset observation, 
where the Galactic component (b) is prominent, with the 
[ apec (cyan) + apec (orange) + phabs x vapec (magenta) ] 
model, using the ARF response, A B (Ei). 

In a strict sense, using A B for the component (b) is not 
correct, which should be A v instead. We made this choice 
due to the limitation of the Xspec vll (§4.1), however, 
it does not matter practically if we only notice the shape 
of the spectrum. The absolute surface brightness of the 
Galactic component was evaluated separately using the 
Xspec fake command and the A v response. We then 
fitted the central region, fixing the shape of the Galactic 
component, but with its normalization scaled so that the 
surface brightness is preserved between the two different 
sky regions. Xspec v12 can handle this situation more 
straightforwardly. 

The released RMF, ae_xi [0-3] .20060213 . rmf, was 
used for the spectral analysis. The ARFs were gener- 
ated by xissimarfgen, and were convolved with the RMFs 
and added for three FI sensors, using the marfrmf and 
addrmf tasks in Ftools. As for the photon statistics 
of the simulation, limit_mode=MiXED was chosen with 
num_photon=100000 and accuracy=0.005. As seen in fig- 
ures 7 (a) and (b), both the observed spectra can be well 
fitted by one temperature plasma emission model for ICM, 
and (apec + apec) model with 1 solar abundance for the 
local Galactic emission. The surface brightness and the 
spectral shape of the Galactic emission is kept constant 
between both regions. This is confirmed by the fact the 
ratio of the Galactic components to the CXB is almost 
equal between figures 7 (a) and (b). Note that the fit ap- 
pears equally good to the BI (XIS1) and FI sensors, which 
are different in the thickness of the contamination. So 
far, it has been confirmed that the temperature for each 
ring derived from the spectral fitting with this method 
is quite consistent with the previous results with XMM- 
Newton and Chandra (Hayakawa et al. 2004; Hayakawa 

10 This value is taken from table 3 of Kushino et al. (2002), for 
the integrated spectrum with source elimination brighter than 
S = 2 X 10~ 13 erg cm" 2 s" 1 (2-10 keV) in the GIS filed of view 
with T = 1.4 (fix) and the nominal NXB level (0%). 
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et al. 2006). See K. Sato et al. in preparation for details 
of the results. 

7. Summary 

• We have developed a Monte-Carlo simulator of the 
Suzaku XRT/XIS system taking into account full 
calibration results. 

• We adopted the ANL platform that provides us 
a flexible and comprehensive environment for the 
Suzaku software production. 

• There is a dedicated task named mkphlist which 
generates a photon file to feed to the simulator. 

• The task xissim reads the photon file, and con- 
ducts the instrumental simulation using the XRT 
ray-tracing library and the RMF of the XIS, and 
generates an event file, which is consistent with 
that for real observation, so that users can analyze 
the simulated data in the same manner as real data. 

• The simulator-based ARF generator is named 
xissimarfgen, which can compute up to 200 ARFs 
corresponding to different accumulation regions by 
a single batch of simulations. 

• The combination of xissim and xissimarfgen enables 
users to analyze spatially extended and spectro- 
scopically complex celestial sources. 

• Since one of the Suzaku's unique features is the 
low and stable particle background, these simu- 
lators are crucial for producing scientific results 
with low signal-to-noise data from extended sources. 

• The latest public version is 2006-08-26, which will 
be included in the next official release of the Suzaku 
Ftools scheduled in late 2006. 
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Table 3. Summary of XIS coordinate column information. 



Column Name 


Min * 


Max T 


Origin 1 


Pixel Size s 


SEGMENT 





3 






rawx/y 





255/1023 




0.024 mm 


actx/y 





1023 




0.024 mm 


detx/y 


1 


1024 


512.5 


0.024 mm 


focx/y 


1 


1536 


768.5 


0.024 mm 


x/y 


I" 


1536 11 


768.5 


0.0002895 deg » 



* TLMINra keywords in the event file. 

I TLMAXn keywords in the event file. 
' TCRPXn keywords in the event file. 
§ TCDLTn keywords in the event file. 

II Default image region. X/Y values can be outside of the region. 
' Angular scale at the center. Outer pixels are slightly different 
due to the tangential projection. 

Table 4. Summary of XIS alignment information 



Item 



Ideal Value 



Focal length 

Optical axis location in DET 
Size of the DET pixel 
Offsets between DET and FOC 
Roll angle between DET and FOC 
Alignment matrix for FOC — ► SKY ' 



4750 mm 
(512.5, 512.5) 
0.024 mm/pixel 
(0.0, 0.0) 
0.0 deg 

3x3 identity matrix 



* Alignment matrix is common to all sensors. 

Appendix 1. Definition of the Coordinates 

The following coordinates are defined to describe event 
locations in the telemetry, on the detector, or on the sky. 

RAW coordinates: Original digitized values in the 
telemetry to identify the pixels of the events. This may 
not reflect physical locations of the pixels on the sensor. 
For example, XIS RAWX (or RAWY) coordinate will have 
values from to 255 (or 1023) on each CCD segment. 
Each of the four XIS sensors has a single CCD chip, and 
a single chip is divided into four segments. 

ACT coordinates: The ACTX/Y values are defined 
to represent actual pixel locations in the CCD chips. 
ACTX/Y will take to 1023 to denote the 1024 x 
1024 pixels in the chip. The XIS RAW to ACT conver- 
sion depends on the observation modes (such as Window 
Options) and will require housekeeping information. The 
XIS ACT coordinate is defined by looking down on the 
sensors, hence the ACTX/Y to DETX/Y conversion needs 
a flip in the Y-direction. 

DET coordinates: Physical positions of the pixels 
within each sensor, XISO-3. Misalignments between the 
sensors are not taken into account. The DETX/Y coordi- 
nate are defined by looking up the sensor, such that the 
spacecraft (S/C) +Y direction becomes the — DETY di- 
rection (the same convention as with ASCA ). The S/C 
Z-axis points in the telescope direction, and +Y direc- 
tion is toward the solar paddle. For XIS, the DETX and 
DETY values take 1 to 1024. 

FOC coordinates: Focal plane coordinate common to 
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Fig. 5. (a) Observed Abell 1060 image combined for the central and offset pointings obtained with XISO in the 1-4 keV energy 
range. The image is smoothed with cr = 16 pixel ~ 17" Gaussian, and the estimated NXB and CXB components are subtracted. The 
exposure time is corrected, but vignetting is not corrected. Directions of DETX/Y axes are indicated in the figure, (b) Comparison 
of the observed and the simulated images (1-4 keV) projected to the DETY axis in the offset pointing. The green line shows the 
simulated distribution by xissim assuming an analytical model (double-/? model) obtained with XMM-Newton, and the kT = 3.4 keV 
vapec model spectrum. The blue and red lines show the estimated CXB and NXB distribution, respectively. The black crosses show 
the observed distribution after subtracting the CXB and NXB components, (c) Same as (b), but for the central observation. 




10 16 



Fig. 6. Plots of the XIS1 ARFs for the Abell 1060 cluster of galaxies calculated with limit_mode=MlXED, num_photon=100000, 
accuracy=0.005, and estepfile=DENSE. (a) The primary extension image in DET coordinate (1024 X 1024) for the central observation 
at the projected radius of r < 2'. (b) Same as (a) but for the offset observation at the projected radius of 17' < r < 27'. (c) The 
SPECRESP columns for the central (black) and offset (red) observations plotted against energy. The 90% confidence range for each 
ARF is indicated by cyan or green lines, respectively, (d) The C0NTAMI_TRANSMIS columns for the central (black) and offset (red) 
observations. 



(a) (b) 




0.5 1 2 5 0.5 1 2 5 

Energy (keV) Energy (keV) 



Fig. 7. Example spectra of the Abell 1060 cluster of galaxies, (a) for the central observation, (b) for the offset observation. In both 
figures, red or black crosses represent the observed spectrum with the XIS1 (BI) or XIS0+XIS2+XIS3 (FI) sensor(s), respectively, 
for the upper panels, and the fit residuals for the lower panels. The estimated CXB + NXB spectrum has been subtracted from 
each observed spectrum, and the estimated CXB spectra for XIS1 are indicated by blue crosses. The spectra are fitted by the 
[ apec (cyan) + apec (orange) + phabs X vapec (magenta) ] model in Xspec 11.3.2t indicated by green line for XIS1 and yellow line 
for the FI sensors. The model components are only plotted for the XIS1 spectrum. 
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Fig. 8. Relations between RAWX/Y, ACTX/Y, DETX/Y among the four XIS sensors. The coordinate are defined looking up 
from the XIS toward the XRT. 



all the sensors. Misalignments between the sensors are 
taken into account so that the FOC images of different 
sensors can be superposed. The origin of the FOC coor- 
dinate corresponds to the XIS nominal position for point- 
ing observations. FOC is calculated from DET by lin- 
ear transformation to represent the instrumental misalign- 
ment, i.e., the offset and the roll angle. 

SKY coordinate: Positions of the events on the sky. 
For each XIS event, the equatorial coordinate of the pixel 
center projected on a tangential plane are given. The 
aberration correction due to parallax (i.e., the revolution 
of the Earth around the Sun) is also considered. 

XRT coordinate: This is given by (XRTX, XRTY) in 
mm on the focal plane, or (6, <ft) corresponding to the 
offset angle (') and the azimuth angle (°) with respect to 
the optical axis of each XRT. The location of the optical 
axis on the DET coordinate is defined so that effective 
area of the XRT is maximized. 

The RAW, ACT, DET, FOC and SKY coordinate are 
written in the Suzaku XIS event files. Relations between 
RAWX/Y, ACTX/Y, DETX/Y among the four XIS sen- 
sors are summarized in figure 8. The DETX/Y pixel sizes 
correspond to the physical pixel size of the XIS CCD, 
while the X/Y pixel size corresponds to the angular scale 
of a single CCD pixel at the reference pixel. To allow rota- 
tion of the image and some shift of the pointing direction 
during the observation, the X/Y range is taken slightly 
bigger than y/2 x 1024. The minimum value, maximum 
value, origin of the coordinate (reference pixel location), 
and pixel size are summarized in table 3. 

There is a file called teldef (namely, telescope definition) 
for each sensor. In the primary header of each teldef file, 
alignment data for the individual sensors (DET— >FOC, 
FOC->SKY, and DET— > XRT) are given. The alignment 
parameters in the teldef file are summarized in table 4. In 
the extensions of the teldef files, sensor-dependent addi- 
tional calibration information may be written. For exam- 
ple, the 1st extension of the XRS teldef file has measured 
positions and sizes of the XRS pixels. 

In this scheme, the conversion from RAW to DET does 
not depend on the misalignments between the sensors. 
Therefore, DETX/Y, as well as RAWX/Y, can be writ- 



ten in the event files without having the calibration infor- 
mation. The DET to FOC conversion requires the sensor 
misalignment data. The conversion from FOC to SKY 
is made using the satellite Z-Y-Z Euler angles (eai, ea,2, 
eas) in the attitude file and the 3x3 alignment matrix 
given in the teldef file. One must be careful because this 
conversion is dependent on the observation date and di- 
rection due to the parallax (aberration) correction. The 
magnitude of the correction is about ±20.5" at maximum. 

All the conversions between these coordinates are sup- 
plied in the form of the C functions in the astetool li- 
brary. These functions make use of the information given 
by the teldef file, and it is strongly recommended to use 
them for the coordinate conversions. They are built on 
the ISAS-made mission-independent library named at- 
Functions, which includes basic routines to handle 3- 
dimensional vectors and rotation matrices. There is also 
a frontend of the coordinate conversions in the Suzaku 
Ftools, named aecoordcalc. 

Appendix 2. Structures &: Parameters 

A. 2.1. mkphlist 

The mkphlist task consists of three ANL modules as 
listed in table 5. The SimASTE_Root (we will omit 
SimASTE_ hereafter in the main text) module is a root 
module for the Suzaku simulators, that handles initial- 
ization of random numbers and common CALDB files. 
The PhotonGcn module generates photons according to 
the parameters set by a user, and caches the photon pa- 
rameters in an internal storage area called BNK (Ozaki et 
al. 2006). The PhotonFitsWrite module retrieves the pho- 
ton data from the BNK and writes the data to the photon 
file. By splitting these functions into dedicated ANL mod- 
ules, it is easier to understand the structure of the task, 
and furthermore we can share the modules among several 
tasks. For example, the Root module is used for all the 
SimASTE tasks, and the PhotonGen modules is shared 
with xissim. 

The parameters of the PhotonGen module (table 5) is 
classified into the following five groups: (1) to determine 
the X-ray flux, photonjlux, flux_emin, flux_emax, and geo- 
metricaLarea; (2) to determine the spectral shape of inci- 
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dent X-rays, spec_mode, qdp_spec_file, and energy; (3) to 
determine the spatial distribution on the sky, image_mode, 
ra, dec, sky_r_min, sky_r_max, fits_image_file; (4) to deter- 
mine the photon arrival time to be equal or random inter- 
val steps, time_mode; (5) to determine how many photons 
are to be generated, limit_mode, nphoton, and exposure. 

A. 2. 2. xissim 

Table 6 summarizes the ANL modules and major pa- 
rameters for xissim. It consist of eight modules, the first 
two modules of which are common to mkphlist. 

In the Root module, the simulation_mode, instrume, 
teldef, and leapfile parameters are added (which are ig- 
nored in mkphlist) when compared with table 5. The sim- 
ulation.mode parameter determines the default mode of 
the simulation, and the two defined modes are discard 
and weight. In the DISCARD mode, each absorbed pho- 
ton is discarded, for example, by absorption in the XRT 
thermal shield. In contrast, the weight of the photon is 
decreased by multiplying the transmission probability of 
the thermal shield in the weight mode. The final value 
of the weight is written to the weight column of the 
output event file. This feature enables efficient simulation 
when most of photons disappear during the simulation, 
however one needs to use care in the handling of the simu- 
lation results. The default simulation_mode is DISCARD for 
xissim, whereas simulation_mode=WEiGHT for xissimarf- 
gen to treat the thermal shield transmission separately 
(§4.2). 

The PhotonGen module enables on-the-fly photon gen- 
eration without input photon files, and is usually deac- 
tivated (enable_photongen=no). The same parameters in 
table 5 are usable in this mode. The PhotonRead mod- 
ule reads up to eight photon files, as well as the GTI file 
and the attitude file, and puts the photon data (ra, DEC, 
photon_time, photon_energy) and the Eulcr angles 
at photon_time into BNK. By mixing multiple photon 
files, it is capable of simulating an observation, e.g. hot 
and widely extended emission from a cluster of galaxies 
with cool emission from the core region. 

The ECStoXRTIN module takes care of the pre-XRT 
component. It retrieves the photon data and the Eulcr 
angles, and converts the photon positions into (8, (f>). The 
parallax (aberration) correction and the cos 9 effect are 
also considered here. XRTsim conducts the ray-tracing 
by calling the xrrt library, and the XRTOUTtoDET com- 
pute the detector position hit by the photon. XISRMFsim 
simulates the XIS using the RMF, and XISevtFitsWrite 
write the final output (table 9) into the event file. 

A. 2. 3. xissimarfgen 

Table 7 summarizes the structure and parameters of 
xissimarfgen. It consists of five ANL modules, and three 
out of which are common to xissim. The two dedi- 
cated modules for xissimarfgen are XISarfPhotonGen and 
XISarfBuild, and they closely cooperate to calculate and 
generate the resultant ARF(s) by driving the XRT part 
of the simulator, XRTsim and XRTOUTtoDET. 

In table 7, parameters of common modules to xissim 



Table 5. Structure and parameters of mkphlist. 

Modulo/Parameter * Description 

SimASTE_Root 

(rand_seed) random number seed 

(rand_skip) random number skip count 

SimASTE_PhotonGen 

photonJIux photon flux (photons cm - s _1 ) 

flux_emin lower energy (keV) for photonJIux 

flux_emax upper energy (keV) for photonJIux 

geometricaLarea XRT geometrical area (cm 2 ) 

spec_mode 0:qdp-spec, l:monochrome 

qdp_specjile qdp spectral file for spec_mode=0 

energy energy (keV) for spec_mode=l 

image_mode 0:FITS-image, l:point, 2:uniform 

ra, dec RA, dec (°) for image_mode=l or 2 

sky_r_min min radius (') for image_mode=2 

sky_r_max max radius (') for image_mode=2 

fitsjmagejile image FITS file for image_mode=l 

time_mode 0:constant, l:Poisson 

limit_mode Oinumber of photon, l:exposure time 

nphoton number of photon for limit_mode=0 

exposure exposure time (s) for limit_mode=l 

SimASTE_PhotonFitsWrite 

outfile output photon file name 

* Parameters in parentheses are hidden parameters. 

are omitted, although the simulation_mode parameter is 
set to weight as mentioned in A. 2. 2. We can categorize 
them as follows: (a) to specify the spatial distribution of 
the celestial target on the sky, source_mode, sourceJmage, 
etc; (b) to specify the accumulation region of the de- 
tected events and corresponding output ARF names, re- 
gion.mode, num_region, regfileiV, detmask, andarffileA; (c) 
to specify the photon statistics at each energy, limit_mode, 
num_photon, and accuracy; (d) to specify the energy step 
to calculate the detection efficiency, rmffile and estepfile; 
(e) to specify the observation date and the satellite Euler 
angles, gtifile, date_obs, attitude, eal, ea2, and ea3; (f) to 
specify other calibration information or simulation modes 
or reference of the SKY coordinate, contamifile, aberration, 
aperture_cosine, pointing, reLalpha, reLdelta, and ref.roll. 
Groups (e) and (f) parameters are similar to xissim. 

Group (a) parameters determine the spatial distribu- 
tion of the target on the sky, and one can specify an 
arbitrary FITS image in the SKY or DET coordinate 
(sou rce_mode=SKYFiTS / detfits) . Pixels with negative 
values are treated as zero in the image. Otherwise, a 
location of a point source can be set in the equato- 
rial coordinate in J2000, or SKY- or DET-coordinate 
(source_mode=j2000/SKYXY/DETXY). In addition, a 
uniform-sky emission with respect to the XRT coordinate 
can be selected (source_mode=UNiFORM). When the FITS 
image or the location of the point source is supplied in the 
DET coordinate, its position on the sky will be affected 
by the wobbling of the spacecraft, hence it is not recom- 
mended to use with the attitude file. 

Note that one must specify skyref to use SKY co- 
ordinates. When source_mode=SKYFiTS, skyref is au- 
tomatically read from the FITS header keywords. As 
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Table 6. Structure and parameters of xissim. Table 7. Structure and parameters of xissimarfgcn. 



Module/Parameter * Description 



SimASTE_Root 

instrume instrument (xis0,xis1,xis2,xis3) 

(simulation.mode) 0:DISCARD, 1:weight 
(rand_seed) random number seed 

(rand_skip) random number skip count 

(teldef) teldef file name 

(leapfile) leap second file 

SimASTE_PhotonGen + 

(enable_photongen) enable on-the-fly photon generation 

SimASTE_PhotonRead 

infileTV input photon file(s) up to N = 8 

(gtifile) name of the GTI file or none 

(date_obs) observation start for gtifile=NONE 

(date_end) observation end for gtifile=NONE 

(attitude) name of the attitude file or NONE 

eal, ea2, ea3 Euler angles for attitude=NONE 

(pointing) pointing type, AUTO or USER 

ref_alpha skyref RA (°) for pointing=uSER 

ref.delta skyref DEC (°) for pointing=uSER 

(reLroll) skyref roll (°) for pointing=uSER 

SimASTE_ECStoXRTIN 

(aperture_cosine) consider aperture decrease by cos# 
(aberration) * enable the aberration correction 

SimASTE_XRTsim 

(shieldfile) XRT thermal shield transmission file 

(mirrorfile) XRT mirror geometry file 

(reflectfile) XRT surface reflectivity file 

(backproffile) XRT backside scatter profile file 

SimASTEJCRTOUTtoDET 

SimASTEJCISRMFsim 

xis_rmffile XIS RMF name 

(aberration) * enable the aberration correction 

(xis_contamifile) XIS contamination file or none 
(xis_efficiency) multiply XIS effciency or not 
(xis_chip_select) discard events fallen outside of CCD 

SimASTEJQSevtFitsWrite 

outfile output event file name 

* Parameters in parentheses arc hidden parameters. 

t See table 5 for rest of parameters when enable_photongen=yes. 

t The aberration parameter is read in two modules. 

long as the WCS (world coordinate system; Greisen 
& Calabretta 2002; Calabretta & Greisen 2002) key- 
words are correctly assigned, one may use an image for 
sourceJmage. When source_mode=SKYXY, things arc a lit- 
tle complicated. If pointing=uSER, the ref.alpha, ref.delta, 
and ref_roll parameters arc utilized for skyref. If point- 
ing=AUTO, which is the default, skyref is read from 
the header keywords of the attitude file, ra_nom and 
DECLNOM, and roll of skyref is always set to 0°, unless 
attitude=NONE. If pointing=AUTO and attitude=NONE, 
skyref is calculated from the specified Euler angles, eal, 
ea2, and ea3, as RA = eal, DEC = 90°-ea2, and roll = 0°. 

Group (b) parameters decide the accumulation re- 
gion^) of the detected events. The num_region parameter 
specify the number of regions to be considered in the ARF 
calculation. Up to 200 regions may be specified in a sin- 
gle batch of simulations. One may specify FITS image(s) 



Module/Parameter Description 

SimASTE_Root T 

SimASTE_XISarfPhotonGen 

(pointing) pointing type, AUTO or USER 

reLalpha skyref RA (°) for pointing=uSER 

reLdelta skyref DEC (°) for pointing=uSER 

(reLroll) skyref roll (°) for pointing=uSER 

sourcejnode skyfits,detfits, 

J2000,SKYXY,DETXY, UNIFORM 

sourceJmage FITS image for source_mode=*FrTS 
source_ra RA (°) for source_mode=j2000 

source_dec DEC (°) for source_mode=,l2000 

sources x (pixel) for source_mode=*XY 

source_y y (pixel) for source_mode=*XY 

source_rmin min 8 (') for source_mode=UNlFORM 
source_rmax max 9 (') for source_mode=UNlFORM 
region_mode skyfits,detfits,Skyreg,detreg 
num_region number of accumulation regions 

regfileiV region file #N, N=l ~ num_region 

arffileiV output ARF #N, N=l ~ num.region 

detmask mask image in DET coord, or none 

limit.mode mixed, num_photon, accuracy 

num_photon number of photons for each energy 
accuracy calculation accuracy for each energy 

gtifile name of the GTI file or none 

date_obs date of observation for gtifile=NONE 

attitude name of the attitude file or none 

eal, ea2, ea3 Euler angles for attitude=NONE 
rmffile RMF to retrieve energy bin 

estepfile E step file or dense, medium, sparse 

(contamifile) XIS contamination file or none 
(aberration) enable the aberration correction 
(aperture_cosine) consider aperture decrease by cos 6 

SimASTE_XRTsim f 

SimASTEJCRTOUTtoDET 

SimASTE-XISarfBuild 

t See table 6 for other parameters. 

or DS9-style region file(s) in SKY- or DET-coordinate 
(region_mode=SKYFiTS / detfits / SKYREG / DETREG) . In 
supplying a FITS image, an unbinned image (1536 x 1536 
for SKY, 1024 x 1024 for DET) is needed to avoid am- 
biguity. The skyref is adopted in the same way when 
source_mode=SKYXY, and the header keywords in the 
FITS image(s) are always ignored. One may optionally set 
the detmask parameter to specify a mask image in DET 
or ACT coordinates, which is automatically judged by the 
CTYPEl and CTYPE2 header keywords. The detmask im- 
age is commonly applied to the all of the specified regions, 
so that this feature is useful in excluding the calibration 
source regions, bad CCD columns, and hot/flickering pix- 
els from the accumulation regions. 

In supplying a FITS image with the regfileA^ parame- 
ter, the pixel values are interpreted as follows. After the 
simulation of each photon, the pixel location on the im- 
age is determined. If the pixel value is zero or negative, 
the photon is discarded as a non-detection. If it is posi- 
tive, the weight of the photon is multiplied by the pixel 
value. Therefore, one should normally supply a binary 



No. 
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Table 8. List of columns in the photon file. 



Table 9. List of columns in the simulated event file. 



Column Name Format* 


Unit 


Description 


PHOTON.TIME lD 


s 


arrival time 


PHOTON _ENERGY IE 


keV 


X-ray energy 


RA IE 


deg 


right ascension of incidence 


DEC lE 


deg 


declination of incidence 



* The column format in FITS convention. 'Id' is a double preci- 
sion floating point, and '1b' is a single precision floating point. 

(0/1) mask image, while a gray-scale image can represent 
non-uniform exposure areas. 

Group (c) parameters specify conditions for the 
photon statistics at each simulation energy. When 
limit_mode=NUM_PHOTON, num.photon count of pho- 
tons are generated at each simulation energy, re- 
gardless of the number of detected photons. When 
limit_mode=ACCURACY, photons are generated until the 
relative error of the detection efficiency becomes less 
than the accuracy parameter for all the accumulation re- 
gions. The relative error is calculated by eq. (4). It may 
take so many photons when the detection efficiency is 
quite low in this mode that it is recommended to specify 
limit_mode=MiXED, in which it generates photons until ei- 
ther of the two conditions is fulfilled. If num_photon=0 or 
accuracy=0 in limit_mode=MiXED, the former or the latter 
is ignored, respectively. 

Group (c) parameters arc utilized to specify the simu- 
lation energy step. The XIS RMF specified by the rmf- 
file parameter is only used to retrieve the energy bin, 
Ei, for the output ARF(s), so that an out-of-date RMF 
may do so far as the energy bin is the same. The es- 
tepfile parameter should be one of the preset keywords 
full/dense/medium/sparse, or point to a file that con- 
tains three decimal numbers, -Emm, E max , Ehm, on each 
line. When estepfile=FULL, the simulation is conducted 
at every RMF energy bin (m = 7900 energy steps), ran- 
domizing each photon energy within a ±1 eV range, which 
can result in a very long computation time. When estep- 
file=DENSE/MEDlUM/sPARSE, the simulation is conducted 
at built-in fixed energies of 2303/157/55 steps, respec- 
tively. These energies are optimized considering the edge 
energies, so that even estepfile=SPARSE can produce an 
acceptable quality ARF for the scientific spectral fitting 
analysis with a featureless and moderate-statistics spec- 
trum. 

Appendix 3. Output File Formats 

A. 3.1. Output of mkphlist 

Table 8 denotes the format of the output photon file 
from mkphlist, which is also the input to xissim. The 
time, energy, and direction of the incident photons are 
contained in the photon_time, photon_energy, and 
(ra, dec) columns, respectivly. This file is basically mis- 
sion independent, except for the GEOMAREA keyword in 
the FITS header, which contains the geometrical area of 
XRT (cm 2 ) specified by the geometricaLarea parameter of 
mkphlist. 



Column Name Format* 


Unit 


Description 


1st extension ' 


EVENTS 






TIME 


Id 


s 


detected time 


PHA 


li 


chan 


pulse height (= pi) 


PI 


li 


chan 


pulse invariant 


STATUS 


li 




status flags (= 0) 


GRADE 


li 




event grade (= 0) 


SEGMENT 


li 




CCD segment id 


RAWX 


li 


pixel 


RAW coordinate x value 


RAWY 


li 


pixel 


HAW coordinate y value 


ACTX 


li 


pixel 


ACT coordinate x value 


ACTY 


li 


pixel 


ACT coordinate y value 


DETX 


li 


pixel 


Db 1 coordinate x value 


DETY 


li 


pixel 


Dbl coordinate y value 


FOCX 


li 


pixel 


rUi^ coordinate x value 


FOCY 


li 


pixel 


IOC coordinate y value 


X 


li 


pixel 


SKY coordinate x value 


Y 


li 


pixel 


SKY coordinate y value 


XISX 


lE 


pixel 


floating point value of detx 


XISY 


IE 


pixel 


floating point value of dety 


WEIGHT 


lE 




weight of the event 


PHOTON.TIME 


Id 


s 


copy of input photon file 


PHOTON_ENERGY 


Id 


keV 


copy of input photon file 


RA 


Id 


deg 


copy of input photon file 


DEC 


Id 


deg 


copy of input photon file 


2ND EXTENSION 


'gti' 






START 


Id 


s 


start of valid time 


STOP 


Id 


s 


stop of valid time 



* The column format in FITS convention. 'Id' is a double preci- 
sion floating point, '1b' is a single precision floating point, 'li' is 
a 16-bit integer. 



Table 10. List of columns in the generated ARF. 



Column Name Format* Unit Description 



Primary extension 'wmap' 
[image] Id count 



image of detected events 



1st extension 'specresp' 

ENERG.LO lE keV 

ENERG-HI lE keV 

SPECRESP lE cm 2 

RESPERR lE cm 2 

RESPRERR lE 

XRT_EFFAREA lE cm 2 

SHIELD_TRANSMIS 1e 

CONTAMI.TRANSMIS 1e 

INDEX lj 

S lE 

T lE 

input 1e count 

DETECT lE count 

weisum 1e count 

RELERR lE 



lower energy bin 
higher energy bin 
computed effective area 
error of specresp 
relative error of specresp 
XRT only effective area 
thermal shield transmission 
contamination transmission 
index of simulated energy 
interpolation coefficient 
interpolation coefficient 
number of input photons 
sum of detected events 
weighted sum of events 
relative error of detect 



* The column format in FITS convention. 'Id' is a double preci- 
sion floating point, Te' is a single precision floating point, Tj' is 
a 32-bit integer. 
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Table 11. List of special header keywords in the output ARF. 



Keyword Name 



GEOMAREA 
TELDEF 
LEAPFILE 
CREATOR 

SOURCE_RATIO_REG 

MASK_RATIO_CCD 

MASK_RATIO_REG 

N_PHOTON 

N_DETECT 

N.WEISUM 

RANDSEED 

RANDSKIP 

RANDNGEN 



Description 

geometrical area of XRT (cm 2 ) 
teldef file name 
leap second file name 
xissimarfgen credit and version 
sourcejmage ratio inside selected region 
detmask ratio in the whole CCD area 
detmask ratio inside selected region 
number of input photons generated 
number of events detected * 
weighted sum of events detected * 
random number seed, = rand_seed 
random number skip count, = rand_skip 
number of random numbers generated 



* Events fallen on a pixel with a positive value for at least one of 
the accumulation regions are treated as "detection" , here. 



A. 3. 2. Output of xissim 

Table 9 denotes the format of the output event file 
from xissim. The output is a standard FITS event file 
with events and GTI extensions, plus information on the 
faked input photons. If the attitude file is given, posi- 
tion on the sky (x and y) is corrected for the wobbling of 
the spacecraft, as well as the parallax (aberration). For 
the calculation of the X and Y columns, the sky reference 
(skyref) of (RA re f , DEC re f, Roll re f) must be specified by 
users as the parameters of ref .alpha, ref .delta, and ref.roll, 
although ref_roll is hidden and usually set to 0. The (x, y) 
positions have a one-to-one correspondence to sky direc- 
tions, however these direction do not necessarily coincide 
with the incident direction of the photon, (ra, dec), due 
to blurring by the PSF of the XRT. 

A. 3. 3. Output of xissimarfgen 

Table 10 summarizes the format of the output ARF gen- 
erated by xissimarfgen. In comparison with the minimum 
set of the ARF (energ_lo, energ_hi, and specresp 
columns), it has several additional columns and the pri- 
mary image extension to record the information of the 
simulation. All of the parameters for xissimarfgen are 
written in the history of the output ARF, and several im- 
portant simulation parameters and results are written to 
the FITS header keywords listed in table 11. 

The energ_lo and energ_hi columns, both in units 
of keV, are copied from the input RMF, and xissimarf- 
gen assumes Ei = (energ_lo + energ_hi) / 2 and AE t = 
energ_hi — energ_lo for the i-th row (i = 1 <~ ra, and 
to = 7900 for the nominal RMF). The specresp column 
holds the final result of the detection efficiency times the 
geometrical area, S A(Ei) 1 in unit of cm 2 . The value of S 
(cm 2 ) assumed in the calculation is written to the FITS 
header keyword of GEOMAREA (table 11). The resperr 
and resprerr columns hold the absolute and relative 
errors, respectively, estimated for the specresp column, 
i.e. resperr = resprerr x specresp. Details are de- 
scribed in § 4.2 for the XRT_EFFAREA, SHIELD_TRANSMIS, 

and CONTAMI_transmis columns. 



The input, detect, and weisum columns hold inter- 
polated values of jV; n , Ad t, arid N w , respectively. The 
weisum represents the weighted sum of events, and is 
usually equal to detect unless one supplies a gray-scale 
region file. The relerr is an interpolated value of the 
relative error calculated by eq. (4) . The coefficients of the 
interpolation, Si and T i; are defined as, 

V(E t ) = s, V(E' mBE J + Ti mU+i)> ( Al ) 

S i = (^Wx, + 1 - E i) I ( E LzX t + l - ^INDExJi ( A2 ) 

t, = (Ei - E' mBSK . ) / (E' mBSKi+1 - ^LJ, (A3) 

where E' mlsx . denotes the simulated energy indexed by 
an integer value of index^, V(Ei) and V(E mlllsx .) de- 
notes a value at the i-th row and the simulation value 
at E mDEXi keV, and Si, Ti, and iNDEXi denote the i-th row 
values of the s, t, and index columns, respectively. As 
one can see easily from these formulae, a simple linear in- 
terpolation is adopted in xissimarfgen. Note that these 
column values in this paragraph do not include the ef- 
fect of the XIS contamination, so that relerr is slightly 
different from resprerr. 

There is a primary extension image written in the out- 
put ARF, too. This image holds the collection of all the 
simulated photons detected in the specified accumulation 
region for the ARF calculation. The coordinates of the 
image are dependent on the region.mode; SKY coordinate 
(1536 x 1536) for region_mode=SKY*, and DET coordi- 
nate (1024 x 1024) for region_mode=DET*. The weight 
value of each photon without contamination is filled in the 
image, so that the image usually contains integer pixel val- 
ues, unless one supplies an gray-scale region file. Pixels 
out of the accumulation region are filled with a value of 
— 1, which is useful in checking that the accumulation re- 
gion is correctly assigned. One can also examine whether 
the celestial target is correctly placed in the FOV. 
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